Intrinsic anomalous Hall effect in nickel: An GGA+U study 
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The electronic structure and intrinsic anomalous Hall conductivity of nickel have been calcu- 
lated based on the generalized gradient approximation (GGA) plus on-site Coulomb interaction 
(GGA-I-U) scheme. The highly accurate all-electron full-potential linearized augmented plane wave 
method is used. It is found that the intrinsic anomalous Hall conductivity io-^y) obtained from the 
GGA+U calculations with U — 1.9 eV and J = 1.2 eV, is in nearly perfect agreement with that 
measured recently at low temperatures while, in contrast, the a^y from the GGA calculations is 
about 100 % larger than the measured one. This indicates that, as for the other spin-orbit interac- 
tion (SOI)-induced phenomena in 3d itinerant magnets such as the orbital magnetic magnetization 
and magnetocrystalline anisotropy, the on-site electron-electron correlation, though moderate only, 
should be taken into account properly in order to get the correct anomalous Hall conductivity. The 
intrinsic a^y and the number of valence electrons (A^e) have also been calculated as a function of 
the Fermi energy (Ef)- A sign change is predicted at Ef = —0.38 eV (A''e — 9.57), and this explain 
qualitatively why the theoretical and experimental a^y values for Fe and Co are positive. It is also 
predicted that fee Ni(i_j;)Co(Fe,Cu)i alloys with x being small, would also have the negative a^y 
with the magnitude being in the range of 500 ~ 1400 f2~^cm~^. The most pronounced effect of 
including the on-site Coulomb interaction is that all the d-dominant bands are lowered in energy 
relative to the Ef by about 0.3 eV, and consequently, the small minority spin X2 hole pocket dis- 
appears. The presence of the small X2 hole pocket in the GGA calculations is attributed to be 
responsible for the large discrepancy in the a^y between theory and experiment. 

PACS numbers: 71. 15. Mb, 71.70.Ej, 72.15.Gd, 75.47.Nb 



I. INTRODUCTION 

Anomalous Hall effect (AHE) refers to the trans- 
verse charge current generation in solids in a ferromag- 
netic phase by the electric field and has received inten- 
sive renewed interest in recent years mainly because of 
its close connection with spin transport phenomena[l|. 
There are several competing mechanisms proposed for 
the AHE. Extrinsic mechanisms of skew scattering [2| and 
side jump[3| refer to the modified impurity scattering 
caused the spin-orbit interaction (SOI). Another mech- 
anism arises from the transverse velocity of the Bloch 
electrons induced by the SOI, discovered by Karplus and 
Luttinger[^], and thus is of intrinsic nature. This in- 
trinsic AHE has recently been reinterpreted in terms of 
the Berry curvature of the occupied Bloch states. p-Q 
Furthermore, recent quantitative first-principles studies 
based on the Berry phase formalism showed that the 
intrinsic AHE is important in various materials. [9| In 
particular, in itinerant ferromagnets such as Fe and Co, 
the intrinsic anomalous Hall conductivity given by first- 
principles calculations^ J^] has been found to agree with 
the experimental anomalous Hall conductivity |ll| - [l3| 
within 30 %, thereby demonstrating the dominance of 
the intrinsic mechanism. Nonetheless, the physical ori- 
gin of the AHE in nickel is still not fully understood. Re- 
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cent first-principles density functional calculations with 
the generalized gradient approximation (GGA) predicted 
a large intrinsic anomalous Hall conductivity of -2203 
i7~^cm~^ in Ni[lO[, which is more than three times 
larger than the corresponding experimental value of - 
646 fl~^cin~^ [1^ at room temperature. In the latest 
experiment 19], the intrinsic Hall conductivity was found 
to be -1100 i7~^cm~^ at low temperatures. Though this 
value is significantly larger than that of the much ear- 
lier experiment [14|, it is still only half of the calculated 
intrinsic Hall conductivity |10l| . 

First-principles GGA calculations have been rather 
successful in describing many physical properties such 
as crystal structure, elastic constant and spin magnetic 
moment, of itinerant ferromagnets Fe, Co and Ni (see, 
e.g., Ref. [la and references therein). However, GGA cal- 
culations fail in describing some relativistic SOTinduced 
phenomena in these itinerant magnets. For example, the 
theoretical values of orbital magnetic moment account 
for only about 50 % of the measured ones in Fe and 
Co jig and the calculated magnetocrystalline anisotropy 
energy of Ni is even wrong in signjlS]. This failure of 
the GGA is generally attributed to its incorrect treat- 
ment of the moderate 3(i electron-electron correlation in 
these systems. Several theoretical methods that go be- 
yond the density functional theory (DFT) with the lo- 
cal density approximation (LDA) or GGA, such as the 
orbital-polarization correction[l^ and LDA/GG A p lus 
on-site Coulomb interaction U (LDA/GGA-FU)[23422 
schemes, have been developed for better description of 



the SOI-induced phenomena in magnetic sohds. In- 
deed, the orbital-polarization correction has been found 
to bring the calculated orbital moments in many itiner- 
ant magnets such as Fe and Co in good agreement with 
experiments [1 71] . Furthermore, it has been demonstrated 
that the correct easy axes and the magnitudes of the 
magnetocrystalline anisotropy energy of Fe and Ni can 
be obtained within the LDA-I-U scheme. [23] 

Therefore, in this work, we perform GGA-l-U calcula- 
tions for nickel to better understand the mechanism of 
its anomalous Hall effect. We use highly accurate all- 
electron full-potential linearized augmented plane wave 
(FLAPW) method(2j]. We find that including on-site 
electron-electron correlation in nickel has significant ef- 
fect on its anomalous Hall effect as well as its electronic 
structure near the Fermi level. In particular, the cal- 
culated anomalous Hall conductivity reduces from -2200 
fi-^cm-i (GGA) to -1066 fi-^cm-^ and the latter value 
is in very good agreement with the measured low temper- 
ature intrinsic Hall conductivity of -1100 J7^^cm^^ from 
latest experiments. [l5| 

This paper is organized as follows. In next section, 
we briefly describe how the intrinsic anomalous Hall con- 
ductivity is calculated within the linear-response Kubo 
formalism as well as the numerical method and compu- 
tational details used in the present work. In Sec. HI, we 
first present the calculated anomalous Hall conductiv- 
ity, magnetic moments and also relativistic band struc- 
ture. We then compare our results with available ex- 
periments and also previous calculations. Finally, we 
make some predictions about the intrinsic anomalous 
Hall conductivity for fee Ni(i_3,)Co(Fe,Cu)2; alloys with 
small Co(Fe,Cu) concentration x, within the rigid band 
approximation. In Sec. IV, we summarize the main con- 
clusions drawn from the present work. 



II. THEORY AND COMPUTATIONAL DETAILS 

The intrinsic anomalous Hall conductivity of a solid 
can be evaluated by using the Kubo formula (2^. The 
intrinsic Hall effect comes from the static limit {uj=0)_oi 
the off-diagonal element of the optical conductivity [2^, 
l26j . Following the procedure for the calculation of the 
intrinsic spin Hall conductivity [2g, we first calculate the 
imaginary part of the off-diagonal element of the optical 
conductivity 



then obtain the real part from the imaginary part by a 
Kramers-Kronig transformation 
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where Vc is the unit cell volume, huj is the photon energy, 
]kn > is the nth Bloch state with crystal momentum k, 
Vx(y) is the velocity operator, and e„'„ = e^.^' — ekn. We 



4',^^) = - 






r^iu) 



(2) 



where P denotes the principle value of the integral. The 
intrinsic anomalous Hall conductivity a^ is the static 
limit of the off-diagonal element of the optical conduc- 
tivity cri j/(w=0). We notice that the anomalous Hall 
conductivity of bcc Fe calculated in this way[^ [231 is in 
quantitative agreement with that calculated directly by 
accounting for the Berry phase correction to the group 
velocity [9,]. 

Since all the intrinsic Hall effects are caused by the 
SOI, first-principles calculations must be based on a rel- 
ativistic band theory. Here the relativistic band structure 
of fee Ni is calculated using the highly accurate FLAPW 
method, as implemented in the WIEN2K codef^SJ. The 
self-consistent electronic structure calculations are based 
on the DFT with the GGA for the exchange correlation 
potential [29|] . To further take d-electron correlation into 
account, we include on-site Coulomb interaction U in 
the GGA-fU approach 20]. The so-called around-mean- 
field (AMF) scheme for double counting correction [22,] is 
adopted here. C/ = 1.9 eV and J = 1.2 eV, which were 
found to give the correct sign and magnitude of magne- 
tocrystalline anisotropy energy for fee Nij23], are used. 
Furthermore, the GGA-l-U calculations using the dou- 
ble counting correction scheme designed for approximate 
self-interaction correction (SIC) for strongly correlated 
systems such as transition metal oxides [21], have also 
been performed. Nevertheless, the theoretical anomalous 
Hall conductivities from the GGA-l-U calculations us- 
ing both double counting correction schemes, are almost 
identical (see Table I below). Since the AMF scheme 
is believed to be more suitable for metallic systems [23], 
here we will concentrate on the results of the AMF cal- 
culations. Moreover, the GGA-l-U calculations with a 
larger U value of 2.5 eV are also performed to see how 
the variation of U may affect the calculated anomalous 
Hall conductivity. 

The experimental lattice constant a = 3.52 AjSO] is 
used here. The muffin-tin sphere radius {Rmt) used is 
2.2 a.u. The wave function, charge density, and poten- 
tial were expanded in terms of the spherical harmon- 
ics inside the muffin-tin spheres, and the cutoff angu- 
lar momentum {Lmax) used is 10, 6, and 6, respectively. 
The wave function outside the muffin-tin spheres was ex- 
panded in terms of the augmented plane waves, and a 
large number of augmented plane waves (about 70 APWs 
per atom, i.e., the maximum size of the crystal momen- 
tum Kmax — ^1 Rmt) wcrc included in the present calcu- 
lations. The improved tetrahedron method is used for the 
Brillouin-zone integration. [3l| To obtain accurate ground 
state charge density as well as spin and orbital magnetic 
moments, a fine 56 x 56 x 56 grid of 185193 /c-points in 
the first Brillouin zone was used. 
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FIG. 1. (Color online) Calculated anomalous Hall conductiv- 
ity a^y as a function of the inverse of the number of fc-points in 
the Brillouin zone (Nk)- The dashed lines are a polynormial 
fit to the calculated values to get the extropolated value of 
a^y at A^fc = oo. For comparison, the theoretical value by 
Wang et aZ. [10| and the experimental value by Ye et oLJ15| 
are also shown as the horizontal dotted and dot-dashed lines, 
respectively. 



III. RESULTS AND DISCUSSION 

Like the calculation of the magnetocrystalline 
anisotropy energy of bulk magnets [l8|, a very fine k- 
point mesh is needed for the anomalous Hall conduc- 
tivity calculation.!^ [13] Therefore, we perform the a^ 
calculations using several extremely fine fc-point meshes 
with the finest /c-point mesh being 106 x 106 x 106. The 
calculated a^ is plotted as a function of the inverse of 
the number (Nk) of fc-points in the first Brillouin zone 
in Fig. 1. The calculated values of a^y are fitted to a 
polynormal to get the converged theoretical a^ (i.e., the 



extropolated value of cr^ at Nk = oo) (see Fig. 1). The 
theoretical ct^ obtained in this way as well as the calcu- 
lated spin {nis) and orbital (nio) magnetic moments are 
listed in Table I. Also listed in Table I are the available 
experimental and previous theoretical <7^ , iris and nio- 

Table I shows that the theoretical a^., from the GGA 
is -2200 17^^cm^^, and is much larger than the experi- 



mental value of -646 fl^^cm"-'^ at room temperature 14 
and also the experimentally derived intrinsic value of - 
1100 ri~-'^cm~-'^ at low temperatures _^15|. Nevertheless, 
the present result is in very good agreement with the pre- 
vious GGA calculations [13, as it should be (see Table I 
and Fig. 1). Interestingly, when the on-site Coulomb in- 
teraction is taken into account via the GGA-f-U scheme, 
the calculated a^y reduces significantly. In particular, 
when U — 1.9 eV and J = 1.2 eV which were found 
to give rise to the correct easy axis and magnitude of 
the magnetocrystalline anisotropy energy [23|, the theo- 
retical a^y becomes -1066 Q~^cui~^. This good agree- 
ment between the present GGA-I-U calculation and the 
low temperature measurements [15] (Table I and Fig. 1) 
indicates that the intrinsic AHE dominates in nickel and 



TABLE 1. Calculated anomalous Hall conductivity a^y 
(n^^cm"^) as well as spin magnetic moment rris (/is/atom) 
and orbital magnetic moment rUs (/is/atom). Superscripts of 
SIC indicate that the values are obtained from the GGA+U 
calculations with the SIC double counting correction scheme 
(see Sec. II). The corresponding experimental values as well 
as previous theoretical a^y are also listed for comparison. 
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that the Sd electron-electron correlation in itinerant mag- 
nets such as nickel has an important effect on the AHE, 
although being only moderate. Further increasing U to 
2.5 eV reduces the calculated a^ slightly (Table I and 
Fig. 1). As expected, the GGA+U calculations increase 
the theoretical spin and orbital magnetic moments, and 
hence enlarge somewhat the discrepancies between the 
calculations and experiments (Table I). Nonetheless, with 
U = 1.9 eV and J = 1.2 eV, the theoretical spin and or- 
bital magnetic moments are still in reasonable agreement 
with the experiments (Table I). 

To help understand how the on-site electron-electron 
correlation affects the electronic band structure and AHE 
in nickel, we plot in Fig. 2 and Fig. 3 the relativistic en- 
ergy bands along the high symmetry lines in the Brillouin 
zone calculated both without and with on-site Coulomb 
interaction U. The relativistic band structure may be 
regarded as the result of a superposition of the corre- 
sponding scalar-relativsitic spin-up and spin-down band 
structures with many accidental band-crossings (degen- 
eracies) lifted by the SOI. In nickel, nevertheless, these 
SOTinduced splittings are generally much smaller than 
the exchange splittings, and thus can be treated as a 
perturbation. [33| Also, including the SOI would lower 
the symmetry of the system. In the present work, the 
magnetization is assumed to be along the [001] direction, 
and the symmetry of the system becomes the tetragonal 
one. The six high symmetry X points which are equiv- 
alent in the nonrelativistic (or scalar-relativistic) case, 
now form two inequivalent groups, namely, two equiva- 
lent ± Z points and four equivalent ± X and ± Y points. 
Therefore, the energy bands near the Z point arc slightly 
different from that near the X and Y points, as shown 
in Fig. 2 and Fig. 3. The relativistic band structure 
and Fermi surface of Ni has been reported by several 
researchers before by using different band structure cal- 
culation methods (see, e.g., Refs. [l3,[la, [23, and [331). In 
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FIG. 2. (Color online) Relativistic GGA and GGA+U band structures. The GGA+U band structure was obtained using 
U = 1.9 eV and J = 1.2 eV. The Fermi level (dotted horizontal line) is at eV. The thick curves denote band 8's which are 
also indicated by the arrowed lines. 



particular, Wang and Callaway analysed in detail the en- 
ergy band characters and Fermi surface sheets. [3j] The 
present GGA relativistic band structure (the red dashed 
curves in Fig. 2 and Fig. 3) is very similar to that 
reported by Wang and Callaway [33'|. For example, in 
both cases, there are five bands (bands 8-12) crossing 
the Fermi level, and there is a small down-spin-dominant 
hole pocket (band 8, noted as X2) centered at the X(Y,Z) 
symmetry point (see Figs. 2-3, and also Fig. 1 in Ref. 
I33I ). However, the small X2 hole pocket was not found in 
the de Haas- van Alphen experiments. [3J) 

Now let us focus on the changes in the relativistic band 
structure caused by including the on-site Coulomb inter- 
action. A pronounced change is that all the d-dominant 
bands are lowered in energy relative to the Fermi level, 
when the on-site Coulomb interaction is taken into ac- 
count in the GGA+U scheme (see Figs. 2-3). In other 
words, the binding energy of the d-dominated valence 
bands is increased by about 0.3 eV. In particular, the 
down-spin band 8 is now pushed completely below the 
Fermi level and the small X2 hole pocket disappears. 
The absence of the X2 hole pocket is consistent with 
the de Haas- van Alphen experiments. 34 [ In Ref. [2j, 
the absence of the X2 hole pocket caused by the on-site 
Coulomb interaction was regarded as the main reason 
that the LDA+U calculations predicted the correct easy 
axis for nickel. As will be shown below, this band 8 near 
the Fermi level calculated without the on-site Coulomb 
interaction, gives rise to a pronounced contribution to 
the anomalous Hall conductivity, and thus the absence 
of the small X2 hole pocket in the GGA+U band struc- 
ture is also the main cause for the significant reduction 
of the calculated anomalous Hall conductivity. 

In Fig. 3, we display the anomalous Hall conductivity 



{(T^y) and the number of valence electrons (A'e) as a func- 
tion of the Fermi level {Ep), together with the relativis- 
tic band structure, from the GGA+U calculations using 
U = 1.9 cV and J = 1.2 cV. The fine /c-point mesh of 
106 X 106 X 106 is used. Clearly, the magnitude of the a^ 
peaks just above the true Fermi level {Ep = eV), with 
a large value of -1420 f^-^cm-i at 0.16 eV (iVg =- 10.3). 
The peak may be related to the fiat band (band 10) along 
the W-L-F line near the L point just above the Fermi 
level (see Fig. 3a). As the Fermi level is further artifi- 
cially raised, the size of the cr^y decreases steadily and 
becomes rather small (within 200 0~^cm~-'^) above 0.75 
eV. When the Ep is artificially lowered, the magnitude 
of the cr^y initially decreases gradually, and then drops 
sharply starting at Ep = —0.35 eV. The resultant shoul- 
der at -0.35 eV can be attributed to the presence of the 
top segment of band 8 near the X (Y,Z) point in this en- 
ergy region (see Fig. 3a). The a^ then changes its sign 
at -0.38 eV (or N^ =~ 9.57). As the Ep is further low- 
ered, the (T^ increases sharply and then peaks at -0.74 
eV (or Ne =- 8.90) with a large value of 2635 r^^icm^i. 
Beyond this point, the cr^ decreases and fiuctuates but 
remains positive as the Ep is further lowered, and it then 
changes its sign again at -2.28 eV (or N^ —^ 4.78). For 
the Ep being below -3.80 eV, the magnitude of the a^ 
is small. Note that the energy bands below -4.0 eV and 
above 0.5 eV are predominantly of 4s4p character and 
thus the effects of the SOI and exchange interaction are 
small. 

Experimentally, the measured anomalous Hall resis- 
tivity pxy is often analysized in terms of two distinctly 
different resistivity (p2;a:)-dependent terms|l|, i.e., pxy = 
apxx + bp^.^ . Since usually p^y -C Pxx , the total anoma- 
lous Hall conductivity 0-3;^ = ~Pxy/{ply+plx) ~ a^xx+b, 
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FIG. 3. (Color online) Relativistic band structure (a), anomalous Hall conductivity (a^y) (b), and number of valence electrons 
(A^e) (c), calculated with (solid blue lines) and without (dashed red lines) on-site Coulomb interaction U. Both a^y and A^e 
were calculated using the fine fc-point mesh of 106 x 106 x 106. The GGA+U results were obtained using [/ = 1.9 eV and 
J = 1.2 eV. The Fermi level (dotted horizontal line) is at eV. In (a), the thick curves denote band 8's which are also indicated 
by the arrowed lines. 



where the first linear cra;a;-dependent term (aa^x) was 
attributed to the extrinsic skew scattering mechanism 
(^fif)0- The skew scattering contribution has been 
found to become dominant in dilute impurity metals at 
low temperatures.!] Indeed, recent ab initio calculations 
for the alloy systems Fci-^Pd^: and Nii_2;Pda; indicated 
that in the small Pd concentration x region, the a^^ 
can be several times larger than the b, and can also dif- 



fer in sign. 35| The second scattering-independent term b 



was usually further separated into the intrinsic contribu- 
tion (T^[4| which can be obtained from band structure 
calculations Ij, as done here for nickel, and also the ex- 
trinsic side jump mechanism (o'|y)[3|. Therefore, when 
one compares the calculated tr^ with the experimental 
scattering-independent term b, one should be aware of the 
possible side jump contribution a^^ even though it has 
been shown that the a^ would dominate in ferromagnet 
metals such as Fe and Co[3,[l0|. Recent theoretical cal- 
culations for the 2D Rashba and 3D Luttinger Hamiltoni- 
ans using a Gaussian disorder model potential suggested 
that the AHE in the (III,Mn)V ferromagnetic semicon- 
ductors at low temperatures could be dominated by the 
gf^.[36| However, more recent ab initio calculations show 
that in fee Fci-^^Pd^r and Nii^jjPdx, the a^y is generally 
two-order of magnitude smaller than both the a^^ and 
a^y, thus being negligible. [35| This is consistent with the 
recent experimental finding of the negligible a^y in bcc 
Fe|13| and fee Ni[i^ using a newly established empirical 
axx-scaling formula for g^„ [13| . 

Let us now return to the calculated a^,, as a function 
of the Ep described above and also presented in Fig. 
3. We find that it can explain qualitatively within the 



rigid band approximation the known AHE experiments 
on 3d transition metal ferromagnets and their alloys. [ll|- 
[l4| For example, the sign of the a^y of both Fe and Co 
were found to be positive. |lll - [l3| Interestingly, the calcu- 
lated a^y at Ne = 8.0 {Ep = -1.02 eV) is 758 l^-^cm-i, 
being in very good agreement with the previous ab ini- 
tio calculations for Fe[9|. Of course, this nearly perfect 
agreement may be accidental because Fe crystallizes in 
the bcc structure and also has a much larger spin mag- 
netic moment, and hence the rigid band model should 
not work very well here. We note that the calculated a^ 
is 2548 r^-^cm-i at N^ = 9.0 {Ep = -0.71 eV), which is 
much larger than the previous theoretical u^ value for 



fee Co [331 and also the experimental tr^ value for poly- 
crystalline hep CoQ. Nickel forms substitutional alloys 
in fee structure with low concentration of Fe, Co and 
Cu. The present calculations (Fig. 3) predict that fee 
Ni(i_2;)Co(Fe,Cu)2; alloys {x being small) would have the 
negative intrinsic tj^ with the magnitude being in the 
range of 500 ~ 1400 il~^cin~^. Of course, when com- 
paring this prediction with the AHE experiments on fee 
Ni(i_3,-)Co(Fe,Cu)a; alloys, one should take into accoint 
the non- negligible skew scattering contribution gf^lSSf. 
Indeed, as mentioned earlier. Fig. 3b shows a sign change 
of crf^ at Ne w 9.570 {Ep = -0.38 eV), whereas the sign 
change of the axy in fee Nii_a;Fea; alloys with x ~ 0.13 
(i.e., Ne ~ 9.74) was experimentally found[3]. This ex- 
perimentally found sign change at the lower Fe concen- 
tration could be attributed to the presence of the positive 
a^^ in fee Nii^^jFe^ alloys. 

As mentioned before, the most pronounced effect of 
including the on-site Coulomb interaction is that all the 



d-dominant bands are lowered in energy relative to the 
Ep by about 0.3 eV. Consequently, the down-spin band 
8 is now pushed completely below the Ep and the small 
X2 hole pocket disappears. Now the top of band 8 is 
located at -0.20 eV on the X(Y,Z) point (Figs. 2-3). To 
see clearly the effect of the on-site Coulomb interaction 
on the cr^ , we also display in Fig. 3 the anomalous Hall 
conductivity (cr^) and the number of valence electrons 
(iVe) as a function of Ep, together with the relativis- 
tic band structure, from the GGA calculations using the 
fine fe-point mesh of 106 x 106 x 106. Strikingly, the 
a^y exhibits a pronounced negative peak right at the Ep 
with the value being about -2178 ri~^cm~-'^ (see Fig. 3b). 
The remarkable difference in the a^., at 0.0 eV from the 
GGA and GGA-I-U calculations can be attributed to the 
presence of the above mentioned X2 hole pocket in the 
GGA calculation (Fig. 3). Another major difference is 
the absence of the clear shoulder at -0.35 eV in the GGA 
a^ spectrum. This is because of the absence of the top 
segment of band 8 near the X (Y,Z) point in this en- 
ergy region in the GGA band structure (see Fig. 3a). 
In fact. Fig. 3 shows that as the Ep is lowered from 
the true Fermi level (0 eV), the GGA a^^ decreases in 
magnitude steeply to zero and changes its sign at -0.30 
eV {Ne = 9.55). Nevertheless, in both cases, the tr^ 
changes its sign near Ne = 9.56. The other pronounced 
difference (by about 700 ri~^cm~^) in the two calculated 
a^y spectra is the height of the positive peak near -0.70 
eV. In the rest of the energy region, the two cr^ spectra 
look rather similar. 



low temperatures fTsj . In contrast, the theoretical a^ 
from the GGA-I-U calculations with U = 1.9 eV and 
J = 1.2 eV is in almost perfect agreement with the mea- 
sured one (Table I). This indicates that, as for other SOI- 
induced magnetic phenomena in 3d itinerant magnets 
such as the orbital magnetic magnetization and mag- 
netocrystalline anisotropy, the on-site electron-electron 
correlation, though moderate only, should be taken into 
account properly in order to get the correct anomalous 
Hall conductivity. The most significant effect of includ- 
ing the on-site Coulomb interaction is that all the d- 
dominant bands are lowered in energy relative to the Ep 
by about 0.3 eV. Consequently, the down-spin band 8 is 
now pushed completely below the Ep and the small X2 
hole pocket disappears. The presence of the small X2 
hole pocket in the LDA and GGA calculations is found 
to be the main reason why the large discrepancy in the 
a^y between theory and experiment exists. The intrinsic 
anomalous Hall conductivity (cr^) has also been calcu- 
lated as a function of the Fermi level Ep (or the number 
of valence electrons Ne) in the GGA-I-U scheme. A sign 
change is predicted at Ep = —0.38 eV (or A^e = 9.57), 
and this explain qualitatively why the theoretical and 
experimental ct^ values for Fe and Co are positive. Fi- 
nally, the present calculations (Fig. 3) indicate that fee 
Ni(i_a;)Co(Fe,Cu)a; alloys {x being small) would also have 
the negative intrinsic ct^ with the magnitude being in the 
range of 500 ~ 1400 f^-^cm-i. 



IV. CONCLUSIONS 
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100 % larger than the intrinsic ct^ recently measured at 
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